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ABSTRACT Properties of cumulant- and combinant ratios are 
studied for multihadron final states composed of Poisson dis- 
tributed clusters. The application of these quantities to "de- 
tect" clusters is discussed. For the scaling laws which hold in 
hierarchical clustering models (void scaling, combinant scaling) 
a generalization is provided. It is shown that testing hierarchi- 
cal models is meaningful only for phase-space volumes not larger 
than the characteristic correlation length introduced by Poisson 
superposition. Violation of the scaling laws due to QCD effects 
is predicted. 
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1. Introduction 

In galaxy clustering studies there is accumulating observational evidence in favour 
of the so-called hierarchical models (see e.g. refs. [1-5] and ref. [6] for a historical 
account). Initiated by the work of Carruthers and Sarcevic [7] and the Tucson 
group [8,9] considerable interest has been devoted to this subject over the last few- 
years also in studying the nature of correlations of multihadron final states [10-13]. 
In the field of multiparticle dynamics one speaks about Linked Pair Approximation 
(LPA) because of the special hierarchy of the higher-order correlation functions. 
Throughout this paper we shall use both naming conventions. 

In the framework of the LPA the two-particle cumulant correlations, defined 
in terms of the inclusive single- and two-particle density correlations p\ and p 2 as 

/C 2 (l,2) = p 2 (l,2)- Pl (l) Pl (2), (1.1) 

provide the building blocks of the higher-order cumulants [7,8]. In eq. (1.1) the ar- 
guments denote coordinates e.g. on the rapidity axis. The cumulants(densities) are 
known also as irreducible(reducible) correlations since the lower-order background 
correlation terms involved by p q are subtracted in JC q . The cumulants measure 
genuine (/-particle correlations, the degree of independence of their arguments. In 
the case of full statistical independence (each argument is independent of all the 
others) the (/-particle density p q factorizes into the product of q single-particle 
densities and the (/-particle cumulant IC q vanishes. In the LPA the normalized 
cumulants 

K q (l,...,q) = lC q (l,...,q)/p 1 (l)...p 1 (<l) (1-2) 

are built up as sums of products of linked two-particle normalized cumulants. For 
example, Ks is composed according to 

k 3 (1, 2, 3) = ^ [« 2 (1, 2)«a(2, 3) + k 2 (2, 3)« 2 (3, 1) 

+ «2(3,1)k 2 (1,2)] 

where the constant A 3 is a free parameter to be determined (by definition, A\ = 
A 2 = 1). In the astrophysical literature the constants A q are known as hierarchical 
amplitudes. The Linked Pair Approximation is in fact a special case of the hier- 
archical models [11]. Some types of topologically distinct graphs allowed in the 
construction of higher-order galaxy correlation functions are absent in the LPA. 
These are the so-called non-snake graphs. In lowest order the two schemes are 
equivalent, leading to eq. (1.3), and they coincide at the level of the integrated 
cumulants. 
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Two years before the discovery of the successfulness of eq. (1.3) for galaxy 
correlation data [2] Mandelbrot arrived at the same functional form by constructing 
a simple model for galaxy clustering [14]. In the model the positions of galaxies are 
stepping points of a Rayleigh-Levy random walk which gives rise to an unbounded 
fractal galaxy distribution. Later it was realized that the distribution of galaxies 
cannot be a pure, unbounded fractal; it contradicts to the fine structure of the 
observed galaxy distribution [6]. The model was modified by introducing many 
Rayleigh-Levy random walk fractals distributed according to a Poisson process [1] . 
The Poisson superposition of bounded fractal clumps well reproduces the visual 
appearance of the galaxy distribution in the sky and it is in agreement with the 
hierarchical structure of higher-order correlation functions [1,6]. 

According to recent work of Giovannini, Lupia and Ugoccioni the hierarchical 
structure of cumulant correlations and the Poisson superposition principle play 
an important role in multiparticle dynamics too. The Torino group studies the 
nature of inside-jet correlations in e + e~ annihilations using Jetset 7.2 Monte 
Carlo events in a wide cm. energy range up to yfs = 1000 GeV [13, 15-17]. One 
of the main results is the identification of jets originated by a quark (antiquark) and 
a gluon. The correlations are in agreement with the hierarchical ansatz in quark- 
and gluon-jets. Moreover, it is found that both types of jets consist of Poisson 
distributed clumps of particles called clans. Within gluon-jets the correlations 
are stronger (the clans have larger particle content) which is attributed to QCD 
effects. In gluon initiated jets the dominant mechanism is the self-interaction of 
gluons whereas quark-jets are controlled by gluon bremsstrahlung emission [13]. 

The Poisson superposition introduces a characteristic correlation length in 
hierarchical models. The main goal of the present paper is the study of the conse- 
quences of this newly appearing length scale. Although correlations of multihadron 
final states will be investigated some of the results may have applications in other 
fields, particularly in galaxy clustering studies. 

2. Basic definitions 

Currently available data for hadrons and galaxies do not allow precise determi- 
nation of the higher-order correlation functions. But testing the validity of hi- 
erarchical models is possible through the integrated cumulants that provide the 
factorial cumulant moments of the underlying count distributions. Moments can 
be determined up to 8th order for galaxy catalogs [4,5] and up to 5th order for 
multihadron final states [8]. In this section we collect some basic definitions and 
formulae concerning count distributions and various moments. For more details, 
see e.g. refs. [10] and [18]. 
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Let us start with the generating function for the multiplicity- or count distri- 
butions P n , conveniently defined by 



g{ Z ) = Y J PnZ n . (2.1) 



n=0 



The connection with the density correlation functions p q is given by the power 
expansion of G(z): 

0(*) = 1 + £ L -^r L $i ( 2 - 2 ) 

where the £ g are the integrated densities providing the factorial moments of P n . 
The integrated cumulants IC q , being the factorial cumulant moments f q of P n , are 
defined by the power expansion of the logarithm of Q{z): 

00 ( 7 _ 1 \Q 

lnG(z) = J2—T-U (2-3) 

q=l 

In the analysis of various count distributions occuring in nature the so-called in- 
finitely divisible distributions play a distinguished role. In multiparticle dynamics 
their importance was emphasized by Giovannini and Van Hove [19,20]. A discrete 
distribution is said to be infinitely divisible if its generating function has the prop- 
erty that for all k > integer \jQ{z) is again the generating function of a certain 
distribution [21]. Q(z) satisfies this property if and only if Q(l) = 1 and 

00 

\ng(z)=\ng(0) + J2C q z q - (2-4) 

9=1 

In eq. (2.4) the C q are the combinants of Gyulassy and Kauffmann [22,23]. The 
combinants should obey C q > and 

00 

J^Cq = -In 0(0) < 00 (2.5) 

g=l 

for infinitely divisible distributions. From eq. (2.5) one sees that the probability 
of detecting no particles at all in a certain phase-space volume, the so-called void 
probability, is Po = Q(0) > for discrete distributions satisfying the conditions of 
infinite divisibility. The generating function takes the form 

G{z)=e* V (f^C q {zi-l)\ (2.6) 
^ 0=1 ' 
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in terms of combinants. The C q can be expressed as combinations of the count 
probability ratios V q = P q /Po, 



hence their name. Eq. (2.7) shows that the knowledge of C q requires only a finite 
number of count probabilities. Furthermore, we need not know the probabilities 
themselves: the combinants follow directly from the unnormalized topological cross 
sections since they involve only ratios of probabilities. It is also seen from eq. (2.7) 
that in general the C q can take negative values as well for q > 2. In this case a 
necessary condition of the infinite divisibility of P n is not satisfied. The combinants 
have some advantageous features, e.g. they share common properties with the 
cumulant moments. For further details see refs. [23,24]. 

3. Cumulant- and combinant ratios in 
Poisson cluster models 

In order to study hierarchical models with Poisson superposition we shall utilize 
the so-called Poisson cluster models. These are closely related to the infinitely 
divisible distributions discussed in the previous section. An advantageous feature 
of these models is the absence of free parameters. The clan production picture 
of hadronization developed by Giovannini and Van Hove [19,20] is a well known 
example of Poisson cluster models. We start with the basic equations. 

Assume that the observed events (the phase-space distribution of multihadron 
final states produced in a certain collision process or the galaxy distribution in the 
sky as a single "event") can be decomposed into Poisson distributed clusters. In 
the Poisson cluster models the generating function of the total event multiplicity 
distribution P n takes the form 



where V(z) and H(z) stand for the generating functions of the Poissonian cluster 
distribution and the arbitrary distribution of particles inside the clusters. G{z) 
is obtained as the convolution of V(z) and 7i(z). Discrete distributions having 
a generating function of the above form are known as compound Poisson dis- 
tributions [21]. Eq. (3.1) is another way of writing Q{z) for infinitely divisible 
distributions, eq. (2.6), with C = — ln^(0) and 




(2.7) 



Q{z) = exp (C(H(z) - 1)) 
= V(H(z)) 



(3.1) 




(3.2) 
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In eq. (3.2) the the (/-particle count probability within a single cluster, p q , is found 
to be [24,25] 

C 

Pq = q r > ( 3 - 3 ) 

reminiscent of the definition of count probabilities P n in terms of the topological 
cross sections a n , P n = a n / J2 n a n- 



Since H(0) = each cluster must contain at least one particle, i.e. p = 
C = 0. Accordingly, Q(0) = V(0) and the probability of detecting no particles, 
Pq, provides also the probability of detecting no clusters in a certain phase-space 
volume. Hence the two basic parameters of the model, the average cluster multi- 
plicity, C, and the average multiplicity within the clusters, q, are found to be [12,26] 

C = -lnP and q = n/C (3.4) 

with n = fx = £i being the total event average multiplicity. It is worth noticing 
that eq. (2.6) can be written also in the form 

oo 

S(z) = IJexp(C ff (z«-l)), (3.5) 

9=1 

i.e. as the product of the generating function of a Poisson distribution of particle 
singlets having mean Cx, a Poisson distribution of particle pairs having mean C2, 
and so on [21]. Thus a formally equivalent derivation of the Poisson cluster models 
is possible in which the total events are composed of Poisson distributed particle 
g-tuples with average multiplicities C q instead of identical clusters or clans. This 
view is closely reminiscent of the Mayer cluster expansion of statistical mechanics. 
Recall that the C q can take negative values as well for q > 2. In this case the 
interpretation of combinants as unnormalized count probabilities within identical 
Poisson clusters or as the average multiplicities of Poisson distributed particle 
(/-tuples loses its meaning. 

On the basis of eq. (3.1) let us consider the relationship between quantities 
corresponding to the observed events and quantities at the level of the individual 
clusters. For the total event factorial cumulants we get [27] 



d q 

= C—H(z) 



z=1 (3.6) 
= C( q 

2=1 



with ( q denoting the factorial moments within a single cluster. As obtained in 
eq. (3.3), the combinants of the total events are related to the count probabilities 
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inside the clusters: 



1 d q 
1 d q 



Ca = -<— n ^G{z) 



(3.7) 

-lid*™ z =o = Cp - 
We see that ratios of quantities derived from In Q{z) are equivalent to ratios of 
quantities derived from H(z) and characterize the individual clusters. This feature 
is well suited to "detect" clusters. When both p- and (/-particle correlations exist 
any of the ratios f q j f p , C q /C p or C q j f p should become independent of changes in 
the size of the phase-space volume if a certain characteristic correlation length is 
exceeded. Cuts in phase-space larger than the typical cluster size correspond to 
cuts in the number of clusters and do not influence single cluster statistics. There- 
fore a flattening is expected in the large-scale behaviour of the above quantities, 
provided that translation invariance holds for extended volumes. Flattening of the 
ratios f q /f\ for wide rapidity intervals has already been studied by Dias de Deus 
using the assumption that multihadron final states consist of a fixed number of 
uncorrelated clusters [28]. In this case the constancy follows from the additivity 
property of the cumulant moments. The large binsize (5y > 2) plateau in the 
behaviour of the f q /f± was indeed observed for NA22 data at = 22 GeV and 
for UA5 data at Js = 546 GeV [28,29]. 



4. Hierarchical models, Poisson 
superposition and scaling 

In this section we consider properties of the cumulant- and combinant ratios again 
but the Poisson superposition principle is combined with the simultaneous valid- 
ity of the Linked Pair Approximation. Assuming translation invariance and the 
linked-pair hierarchy of the higher-order cumulant correlation functions the strip- 
integration [8] of K q produces normalized factorial cumulant moments obeying the 
recurrence relation 

K q = kq= A q K\~\ (4.1) 
Jl 

The hierarchical amplitudes A q should be independent of collision energy and 
phase-space volume. Eq. (4.1) is valid also for hierarchical models that allow not 
only snake-type graphs in the construction of higher-order cumulant correlation 
functions. The integrated amplitudes in eq. (4.1) may deviate from the local ones 
appearing in the expressions for K q . In addition to the direct test of eq. (4.1) it 
is possible to check the validity of the LPA through various scaling laws. For the 
ratio of two factorial cumulants eq. (4.1) yields 

4 = ^(f2/fi) q - p - (4-2) 
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Expressing the combinants in terms of factorial cumulant moments [24,25,30], 

r=q 

one sees that the cumulant- and combinant ratios considered in the previous sec- 
tion obey a universal behaviour in hierarchical models. Any of these quantities 
should depend on collision energy and phase-space volume not arbitrarily but only 
through the combination /2//1. Beside eq. (4.2) we have e.g. 



(4.4) 



= Xqpih/fl), 



the combinant-to-cumulant ratios chosen at fixed q, p from different energies and 
volumes scale to a universal curve when plotted against fil f\- In eq. (4.4) the con- 
tribution of an infinite number of factorial cumulant moments is confined into the 
scaling variable fil f\- The Xqp provide the generalization of the scaling functions 
introduced in ref. [25] which are the p = 1 special cases of eq. (4.4). 

The generalization of the famous void scaling function [31,32] can be obtained 
in a similar manner. The probability of detecting no particles in a certain phase- 
space volume obeys the scaling law 

- m -Po \^ (- 1 ) r ~ 1 A r , 
= xopihlh) 

as is seen from eq. (2.3) by substituting z = 0. The zero subscript of x m eq. (4-5) 
refers to the property that these scaling functions involve only the void probability. 
Through eq. (2.5) we get xop = S^li Xqp- Similarly to the void scaling functions, 
the Xqp constructed from the combinants characterize the distribution of empty 
regions in phase-space. The C q can be interpreted as minus the logarithm of the 
probability that a certain phase-space volume is empty of particle g-tuples, see 
eq. (3.5). The advantage of Xqp over Xop lies in the fact that the contribution of 
low-order factorial cumulants can be excluded in a systematic manner by increasing 
q (observe that fi and fi involved by Po carry no information on the validity of 
the LPA). 



Since the relationship between count probabilities and factorial moments is 
the same as the relationship between combinants and factorial cumulants [23,30] 
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a scaling law analogous to eq. (4.4) holds for the count probabilities P n if the 
recurrence relation 

F q =U = VT 1 (4-6) 

SI 

with constant coefficients A q is satisfied for the normalized factorial moments: 



(4.7) 



with n, p > 1 since £o = 1- That is to say, if one plots the ratios P n /Cp chosen 
at fixed n, p against the combination £2/^1 the validity of eq. (4.6) results in a 
universal curve, rj np (^2/^i), instead of many different behaviours corresponding to 
different collision energies and volume sizes. Count probability ratios not involving 
P constitute another set of scaling functions. The analogue of void scaling for 
eq. (4.6) is provided by r) 0p = (1 - Pq)/£ p = Y.n=i Vn P - 

Eq. (4.6) characterizes monofractal density fluctuations. They are expected 
to occur e.g. in second-order phase transitions. Hence the scaling law of the 
multiplicity distributions expressed by eq. (4.7) can in principle signal the forma- 
tion of quark-gluon plasma if QCD undergoes a second-order transition [30]. It 
should be emphasized that eq. (4.6) holds for any field theory with dimensionless 
coupling [33]. For example, in QCD with fixed coupling constant a s the density 
fluctuations are concentrated into randomly distributed monofractal patches of 
phase-space [34]. With running coupling a s multifractal fluctuations appear and 
the exponent on the rhs. of eq. (4.6) gets multiplied by a (/-dependent factor. QCD 
effects violate the above scaling rule of the count probabilities P n . 

We turn our attention to the validity of the LPA and the resulting scaling 
laws if the observed events are composed of Poisson distributed clusters. First of 
all, from eq. (3.6) we get T q = C q ~ x K q and 



with T q denoting the normalized factorial moments within a single cluster, Cg/Ci- 
One sees that the hierarchical amplitudes A q are equivalent to single cluster statis- 
tics. Consequently, they should be independent of changes in the size of the phase- 
space volume if the typical cluster size is exceeded, regardless of the validity of 
the LPA. For the scaling laws we have similar equivalence relations. According 
to eqs. (3.6) and (3.7) the scaling functions of eq. (4.4) at the level of the ob- 
served events are equivalent to the scaling functions of eq. (4.7) at the level of the 
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individual clusters: 

XqyUll fl) I total event = Vqp(&/£l) | gingle cluster - (4-9) 

Moreover, the void scaling functions of eq. (4.5) are equivalent to l/( q , the recip- 
rocal of the unnormalized factorial moments inside the clusters. Thus for multi- 
hadron final states composed of Poisson distributed clusters testing eq. (4.1) and 
the scaling laws expressed by eqs. (4.4) and (4.5) is meaningful only for phase-space 
volumes not exceeding the characteristic correlation length introduced by Poisson 
superposition. Cuts in phase-space larger than the typical cluster size are equiva- 
lent to cuts in the number of clusters and leave single cluster statistics unchanged. 
We mention that the cluster size may vary with energy, e.g. in hh collisions the 
threshold binsize in rapidity at which the constancy of the ratios f q / fi appears is 
increasing from 5y ~ 2 to 5y ~ 3 between = 22 and 546 GeV [28,29]. In the 
X-scaling analyses already performed [11-13] the chosen cuts in rapidity do not 
exceed the above values. 

Closing this section let us reconsider the promising results of the Torino group 
mentioned in the Introduction. The finding that quark- and gluon-jets consist of 
Poisson distributed groups of particles, the clans, and the confirmation of hier- 
archical inside-jet correlations provide a direct connection to our results. The 
hierarchical pattern of cumulants in quark- and gluon-jets was deduced from the 
validity of the void scaling law eq. (4.5) for p = 1 [13]. Due to the Poisson su- 
perposition this scaling rule should follow from the validity of eq. (4.6) at the 
level of the individual clans. Since clan properties inside gluon originated jets are 
controlled by the self-interaction of gluons, deviations are expected from the void- 
and combinant scaling functions in gluon-jets. For running QCD coupling con- 
stant ct s eq. (4.6) is no longer valid [33,34] and this violates the inside-clan scaling 
rules of count probabilities and unnormalized factorial moments being equivalent 
to the inside-jet combinant- and void scaling laws. In the yfs = 1000 GeV cm. 
energy range investigated by the Torino group the self-interaction of gluons may 
cause observable scaling violation for gluon initiated jets. Deviations from the 
void scaling function %oi axe indeed observed [13] and they are attributed to the 
absence of translation invariance. The separation of the two sources of scaling 
violation would be of particular importance. According to preliminary results of 
De Wolf [35] the breakdown of scaling can be seen also in the behaviour of the 
higher-order Xqi combinant-to-cumulant ratios for Jetset Monte Carlo events. 

5. Summary 

One of the common problems in the study of the very large and very small length- 
scale phenomena is the formation of structures — the texture of multihadron final 
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states in phase-space and the distribution of galaxies in the sky. The assumption of 
randomly superimposed objects (clusters, clans) with cumulant correlation func- 
tions obeying hierarchical structure at the level of the total matter distribution 
plays a distinguished role in multiparticle physics and particularly in galaxy clus- 
tering studies. We have investigated the properties of cumulant- and combinant 
ratios for hierarchical models with Poisson superposition. According to our results 
the behaviour of the ratios f q /f p , C q /C p and C q /f p among others considerably 
simplifies with the two assumptions: 

i) From eqs. (3.6) and (3.7) one obtains that for Poisson distributed clusters 
(clans) the above quantities are equivalent to single cluster statistics. Conse- 
quently, they should be independent of changes in the size of the phase-space 
volume if the characteristic correlation length introduced by Poisson super- 
position is exceeded. Cuts in phase-space larger than the typical cluster size 
correspond to cuts in the number of clusters and do not influence single clus- 
ter statistics. Therefore large-scale constancy of the cumulant- and combinant 
ratios is expected at a fixed energy. 

ii) The consequence of the LPA-relation, eq. (4.1), is that the cumulant- and 
combinant ratios should depend on collision energy and phase-space volume 
not arbitrarily but only through the momentum combination f 2 / fi . One ar- 
rives e.g. at the nontrivial scaling laws expressed by eqs. (4.4) and (4.5). 
These are generalizations of the scaling laws introduced in refs. [31] and [25]. 
According to i) for a set of Poisson distributed clusters all of the hierarchical 
statistics (A q , \ qv , xo P ) are equivalent to statistics characterizing the indi- 
vidual clusters. Therefore testing the validity of eq. (4.1) and the resulting 
scaling laws is meaningful only for volume sizes not larger than the typical 
size of the clusters. 

By eq. (4.8) we are led to the conclusion that in hierarchical models with Poisson 
superposition the linked-pair hierarchy of higher-order correlation functions man- 
ifests itself primarily at the level of the individual clusters. Within a single cluster 
the density correlation functions and factorial moments (the reducible statistics) 
obey the hierarchical structure. The cumulant correlation functions and factorial 
cumulant moments (the irreducible statistics) exhibit the hierarchical pattern at 
the level of the total events as the consequence of Poisson superposition. Eq. (4.1) 
and the resulting combinant- and void scaling laws hold for any field theory with 
dimensionless coupling if randomly superimposed clusters are present, e.g. in QCD 
with fixed coupling constant a s . Restoration of the running of a s causes violation 
of the scaling laws. This can be most significant in e + e~ annihilations for gluon 
initiated jets controlled by gluon self-interaction. 
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We close with a list of some measurements that may shed more light on the 
validity of hierarchical models and Poisson superposition: 

• The Poisson cluster models discussed in the present paper can be tested e.g. 
by measuring the C q . Observation of negative combinants signals that the 
underlying count distributions are not infinitely divisible. In this case the 
Poisson cluster models are inappropriate to represent the data. 

• The presence of clusters can be revealed by correlation length measurements. 
They consist of determining the threshold volume size at which the constancy 
of the cumulant- and combinant ratios appears. The threshold provides an 
estimation for the typical size (e.g. rapidity extent) of the Poisson distributed 
clusters. 

• Measurement of the higher-order combinant scaling functions, e.g. the Xqi 
for q > 3, is important to confirm eq. (4.1). These have the advantage over 
the void scaling function xoi that the contribution of fi and is excluded. 
The first two factorial cumulants have considerable influence on the shape of 
Xoi but they carry no information on the validity of eq. (4.1). 

• It would be interesting to see how clear is the violation of ry-scaling if x-scaling 
holds valid. This provides information on the selective power of scaling 
laws that have different origin. The simplest possibility is the comparison 
of — lnP //i and (1 — -Po)/£i as functions of /2//1 and £2/^1 respectively. 

Lack of translation invariance may considerably influence the proposed measure- 
ments. Estimation of its significance is important to gain reliable information on 
higher-order correlations obeying hierarchical structure with the possible presence 
of clusters distributed according to a Poisson process. On the basis of the experi- 
mental data that will be available at LEP200 and TEVATRON, at twice as large 
cm. energies as present, we have the chance of answering the open questions in 
the near future. 
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